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INTRODUCTION 

In  recent  years  a number  of  thermal  mapping  techniques  have 
evolved  into  accepted  diagnostic  tools  to  study  problems  in  aero- 
dynamic heating.  This  report  deals  with  the  so-called  phase-change 
paint  technique,  which  employs  temperature-sensitive  coatings  which 
change  phase,  i.e.,  melt,  from  an  opaque  crystalline  solid  to  a clear 
liquid  at  a specified  temperature.  The  method  was  pioneered  by 
Jones  and  Hunt(l)  at  NASA,  Langley  Research  Center,  in  the  mid- 
sixties. Since  then,  the  technique  has  been  used  at  most  major  wind- 
tunnel  facilities.  (See  Refs.  (2)  and  (3),  for  example.) 

Initially,  the  phase-change  paint  method  was  recognized  as  one 
which  provided  only  a qualitative  look  at  aeroheating  patterns  on 
complex  geometries.  Two  strong  points  of  the  method  were  (1)  that 
the  models  were  generally  easy  to  fabricate,  since  no  extensive 
instrumentation  was  required,  and  (2)  that  these  models  were 
essentially  fully  "instrumented"  for  situations  where  the  interfer- 
ence heating  patterns  were  of  unknown  location  and  extent.  The 
phase-change  paint  technique  lias  since  developed  to  the  point  that 
one  may  expect  to  obtain  reliable  quantitative  data  in  good  agreement 
with  heat-transfer  results  from  methods  employing  thermocouples  or 
heat -flux  gages,  for  example. 

Almost  every  user  of  the  phase-change  paint  technique  generally 
has  seen  fit  to  document  his  results  and  to  comment  on  certain 
aspects  of  the  method.  In  this  respect,  this  report  is  similar  to 
others  in  the  literature.  It  documents  the  use  of  the  phase-change 
paint  technique  at  the  Naval  Surface  Weapons  Center,  White  Oak 
Laboratory,  to  study  aerodynamic  interference  heating  problems  on 
fin-body  configurations.  It  addresses  various  features  of  the 
paint  method  and  presents  some  representative  test  results.  A 
subsequent  technical  report  documents  in  more  detail  the  results  of 
a fin-cone  interference  heating  investigation (^) . 


Jones,  R.  A.,  and  Hunt,  J.  L. , "Use  of  Fusible  Temperature  Indi- 
cators for  Obtaining  Quantitative  Aerodynamic  Heat-Transfer 
Data,"  NASA  TR  R-230,  Feb  1966 

(2) 

'Martindale,  W.  R. , "Interference  Heating  Measurements  on  a Hyper- 
sonic Cruise  Vehicle  Wing  Using  the  Phase-Change  Paint 
Technique,"  AEDC-TR-70-78 , Apr  1970 

^^^Patterson,  J.  L.,  "Heat  Transfer  Testing  in  the  AFFDL  High 

Temperature  Facility  Using  the  Phase  Change  Coating  Technique," 
AFFDL  FXG  TM  70-12,  Aug  1970 

Gillerlain,  J.  D. , Jr.,  "Experimental  Investigation  of  a Fin-Cone 
Interference  Flow  Field  at  Mach  5,"  NSWC/WOL/TR  75-63 
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This  report  elaborates  on  certain  points  of  both  the  underlying 
theory  and  the  experimental  technique  in  an  effort  to  provide  an 
overview  of  the  phase-change  paint  technique.  As  contributions  to 
the  theory  behind  the  method,  analytic  relationships  have  been 
worked  out  which  indicate  what  relative  uncertainty  in  the  heat- 
transfer  coefficient  (as  determined  from  reduction  of  the  phase- 
change  paint  data)  is  produced  by  uncertainties  in  the  various 
input  parameters  relevant  to  the  problem.  This  information  should 
provide  future  users  of  the  method  with  helpful  guidelines  in 
designing  experiments. 


BACKGROUND 

Two  common  types  of  temperature-sensitive  coatings  are  available 
commercially.  Each  indicates  surface  temperature  by  a different 
observable  mechanism.  One  type  is  a so-called  color-change  coating 
(known  as  Detectotemp)  which  exhibits  a specified  color  change  or 
series  cf  color  changes  over  a specified  temperature  range.  The 
second  type  is  a phase-change  coating  (called  Tempilaq)  which  melts 
irreversibly  from  an  opaque  crystalline  solid  to  a clear  liquid  at 
a single  specified  temperature.  A bench-test  program  on  these  two 
types  of  coatings  was  carried  out.  Based  on  results  of  these  lests, 
and  on  results  reported  by  other  researchers,  the  phase-change 
paints  were  selected  over  the  color-change  paints  as  being  better 
suited  for  use  in  a short-duration  high-speed  v;ind-tunnel  facility. 

A series  of  calibration  checks  involving  both  types  of  coatings 
further  reinforced  the  choice  of  Tempilaq  for  our  purposes.  These 
calibration  checks  are  included  as  Appendix  A.  Briefly,  it  was 
found  that  Detectotemp  color-change  coatings  have  a built-in  time- 
temperature  history  associated  with  their  specified  color  changes. 

If  the  color  change  occurs  in  a very  short  time  period  (less  than 
30  minutes) , as  is  usually  the  case  in  high-speed  wind-tunnel  tests, 
then  the  temperature  at  the  instant  of  the  color  change  is  higher 
than  the  rated  temperature.  The  manufacturer  of  Detectotemp  provides 
a brochure (5)  which  contains  calibration  charts  of  color-change 
temperature  versus  heating  time.  In  addition,  these  coatings  may  be 
subject  to  extraneous  color  changes  as  the  result  of  loss  of 
evaporative  salts  when  a vacuum  is  drawn  on  the  wind-tunnel  test  cell. 

The  phase-change  paints  are  considered  to  be  more  reliable 
temperature  indicators  for  use  in  short-duration  (with  a lower  limit 
of  test  times  greater  than  about  two  seconds  for  reasons  to  be  seen 
later)  high-speed  wind-tunnel  tests.  They  are  not  known  to  vary 
with  either  pressure  or  heating  rate  in  melting  at  their  rated 
temperature (1) . Moreover,  they  do  not  involve  the  subjective  element 
of  color  discrimination,  the  melt  process  being  one  of  an  opaque 
solid  turning  to  a clear  liquid. 

Tsl 

Thermochrom  and  Detectotemp  Information  Brochure,  published  by 

H.  V.  Hardman  Co.,  Inc.,  Belleville,  New  Jersey 
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THEORY 


Reduction  of  the  phese-change  paint  heat- transfer  data  is 
based  on  solving  the  transient,  one-dimensional  heat-conduction 
equation  given  by 


3T 

3t 


(1) 


where  T is  temperature,  t is  time,  x is  distance  normal  from  the 

surface  into  the  model,  and  k,  p,  and  C are  the  thermal  conductivity, 

P 

the  density  and  the  specific  heat  of  the  model  material,  respectively. 
The  mathematical  boundary  conditions  and  assumptions  necessary  to 
solve  Equation  (1)  may  be  stated  in  the  following  physical  terms: 


1.  The  model  is  initially  at  a uniform  temperature, 


T(x,0)  = T^  . (2) 

2.  The  depth  of  heat  penetration  into  the  model  is  small 
compared  to  the  wall  thickness  and  surface  radius  of  curvature, 

T(«,t)  = T^  . (3) 

(This  is  the  semi-infinite  slab  approximation.) 

3.  The  local  heat-transfer  coefficient,  h,  is  constant  with 
time,  and  the  heat  flux  leaving  the  gas  at  the  surface  equals  the 
heat  flux  into  the  model, 

k ■■  = h[T^^  - T(0,t)]  , (4) 

where  T^^  is  the  adiabatic  wall  temperature. 

4.  The  model  properties  do  not  vary  significantly  with 

temperature,  so  that  the  thermal  diffusivity,  a = k/pC  , is 
constant  with  time.  ^ 

5.  The  phase-change  paint  coating  is  very  thin,  so  that  it  and 
the  model  surface  are  assumed  to  be  at  the  same  temperature  at  the 
same  time, 

T(0,t)  = , (5) 

where  T^  is  the  rated  temperature  at  which  the  coating  melts. 
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Based  on  the  above  assumptions,  the  solution  (Ref.  (6))  of 
Equation  (1)  may  be  expressed  as 


T = 1 - e^  erfc  3 


(6) 


where  T and  3 are  dimensionless  parameters  defined  to  be 


T - 


T - T. 
m 1 

T - T. 
aw  1 


(7) 


and 


3 = 


/t 


(8) 


and  erfc  3 is  the  complementary  error  function  of  3.  Equation  (6)  is 
shown  graphically  in  Figure  1.  Additional  details  of  the  solution 
of  Equation  (1)  are  included  in  Appendix  B. 


EXPERIMENTAL  TECHNIQUE 


CHOICE  OF  MODEL  MATERIAL 

From  inspection  of  Equations  (6)  and  (8),  it  is  apparent  that 
the  data*  reduction  is  very  sensitive  to  the  thermophysical 
properties  of  the  model.  Teflon  has  proved  to  be  a suitable  model 
material  for  a number  of  reasons: 


1.  It  is  homogeneous  with  known  thermal  properties.  The 
generally  accepted  thermophysical  properties  of  Teflon  are  listed 
in  Table  1. 


2.  It  can  withstand  relatively  high  temperatures. 

3.  Teflon  has  a low  thermal  diffusivity.  Lateral  conduction 
effects  are  minimized.  Thermal  diffusion  times,  t^,  are  large, 

thus  fermitting  adequate  test  times,  the  criterion  being  that  the 
test  time  (the  time  required  for  the  coating  to  melt)  must  be  less 
than  the  thermal  diffusion  time  in  order  that  the  model  behave  like 
a semi-infinite  slab.  The  restriction  on  the  thermal  diffusion  time, 
that  is,  the  time  for  a heat  pulse  to  reach  a given  depth  in  a 
material,  in  order  for  the  semi-infinite  slab  approximation  to  hold 
is  given  by  the  dimensionless  parameter 

^ < 0.2  (9) 


Carslaw,  H.  S., 
Second  edition. 


and  Jaeger,  J.  C.,  Conduction  of  Heat  in  Solids, 
Oxford  Univ.  Press,  Inc.,  1959,  p.  71 
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where  a is  the  thermal  diffusivity  and  S-  is  a characteristic  length, 
which  is  typically  the  model  wall  thickness.  Under  the  restriction 
of  Equation  (9)  the  solutions  for  the  finite  slab  and  semi-infinite 
slab  do  not  differ  by  more  than  10  percent  (Ref.  (1)).  Recently,  a 
number  of  reports  have  provided  information  on  how  to  make  corrections 
for  thin-slab,  i.e.,  finite  slab,  conditions , ^ » 9)  when  such 
corrections  are  deemed  necessary.  The  characteristic  length  in 
Equation  (9)  may  well  be  a radius  of  curvature,  R . Based  on 

solutions  in  the  literature  for  constant  heat  flux  to  a surface, 
the  difference  in  heat-transfer  coefficient  for  a slab  and  cylinder 
does  not  exceed  10  percent  when 


Q.12 


(10) 


thus,  for  this  limiting  condition,  no  curvature  correction  is  needed. 
As  pointed  out  in  Reference  (7) , Equation  (10)  is  only  raiely  a 
more  restrictive  condition  than  Equation  (9) , except  in  a leading- 
edge  region. 

4.  Teflon  IS  sufficiently  strong  to  withstand  loading  from 
rapid  injection  into  the  airstream  (by  means  of  a hydraulic  ram,  for 
example) . 


5.  Teflon  is  impervious  to  the  paint  coating  and  thinner  used 
in  the  tests. 

The  above  factors  support  the  use  of  Teflon  as  a model  material 
in  phase-change  paint  heat-transfer  tests.  The  conclusion  of 
Schultz  in  Reference  (11)  is  considered  to  be  overly  harsh  in 
excluding  Teflon  as  a suitable  material  for  this  technique.  The 


(7) 

Maise,  G.,  and  Rossi,  M.  J. , "Lateral  Conduction  Effects  on  Heat- 
Transfer  Data  Obtained  with  the  Phase-Change  Paint  Technique," 
NASA-CR-2435,  Aug  1974 

/ O \ 

' 'Hunt,  J.  L.,  and  Pitts,  J.  I.,  "Thin  Wing  Corrections  of  Phase- 
Change  Heat-Transfer  Data,"  J.  Spacecraft,  Vol.  8,  1971,  pp. 
1228-1230 

(9) 

Hunt,  J.  L. , Pitts,  J.  I.,  and  Richie,  C.  B. , "Application  of 
Phase-Change  Technique  to  Thin-Sections  with  Heating  on  Both 
Surfaces,"  NASA  TN  D-7193,  Aug  1973 

^^^^Eckert,  E.  R.  G.,  and  Drake,  R.  M. , Heat  and  Mass  Transfer, 

Second  edition,  McGraw-Hill,  Nev;  York,  1959 

^^^^Schultz,  H.  D.,  "Experimental  and  Analytical  Investigation  ff 
Temperature  Sensibive  Paints,"  AFFDL-TR-72-52 , Jun  1972 
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apparent  problem  of  a fluctuation  in  around  90°F  may  be  avoided 
by  raising  the  initial  temperature,  T^,  of  the  model. 


MODEL  CONFIGURATION  FOR  FIN-BODY  INTERFERENCE  HEATING  TESTS 

A fin-body  model  was  constructed  using  an  existing  five-degree 
half-angle  cone  model  made  of  white  Teflon.  A conical  extension  of 
the  same  cone  angle  was  fabricated  from  dark  gray  Teflon.  This 
extension  was  outfitted  with  two  cylindrically  blunted  fins  180° 
apart,  one  unswept  and  one  swept  60°  with  respect  to  the  cone  surface 
normal.  The  fins  were  also  made  of  dark  gray  Teflon.  The  cone  had 
a stainless  steel  nosetip  to  assure  a sharp  tip  and  the  fins  had 
stainless  steel  shaft  inserts  fou  strength.  A photograph  of  the 
model  is  shown  in  Figure  2,  accompanied  by  a schematic  diagram  of 
the  model  in  Figure  3.  The  extension  is  shown  schematically  in 
Figure  4.  Dark  gray  Teflon  was  used  for  the  fin-cone  extension 
because  it  provided  better  color  contrast  between  the  model  surface 
and  the  phase-change  coatings,  most  of  which  dry  to  a light  opaque 
color.  The  extension  was  instrumented  with  four  embedded  thermo- 
couples, as  indicated  in  Figure  4.  These  served  two  purposes. 

First,  they  indicated  the  initial  temperature  of  the  model,  which  is 
a necessary  input  parameter  in  the  data-reduction  scheme  (see 
Eq.  (7)).  Secondly,  they  recorded  any  temperature  rise  at  their 
embedded  depth  in  the  model  resulting  from  thermal  diffusion  of 
surface  heat  transfer,  indicative  that  the  model  was  no  longer 
behaving  thermally  as  a semi-infinite  slab. 


PHASE-CHANGE  COATINGS 

Tempilaq  coatings  with  melt  temperatures  of  163°F,  213°F  and 
263°F  were  selected  for  use  in  most  of  these  tests  based  on 
estimates  of  aeroheating  temperatures.  The  paints  are  available 
over  the  range  100°F  to  2500°F  in  varying  increments  of  three 
degrees  at  the  lower  temperatures  and  50  degrees  at  the  higher  ones. 
The  color  contrast  with  the  model  and  the  crystalline  nature  of  the 
dried  coatings  were  also  factors  in  these  choices. 

The  selected  paints  were  thinned  (using  a special  thinner  manu- 
factured by  Tempilaq)  and  applied  to  the  model  by  means  of  an  air- 
brush. Careful  use  of  an  airbrush  permits  one  to  obtain  a thin, 
fairly  uniform  coating.  Only  a very  thin  coating  (less  than  1 mil) 
is  needed,  the  assumption  being  that  it  and  the  surface  are  at  the 
same  temperature  at  the  stme  time  (Eq.  (5) ) . Too  heavy  a coating 
may  result  in  running  of  the  coating  as  it  melts.  In  addition,  too 
thick  a coating  could  introduce  error  due  to  latent  heat  of  melting 
which  is  not  accounted  for  in  the  theory. 


DATA  ACQUISITION 

The  raw  data  from  these  tests  consisted  of  color  16mm  motion 
pictures  of  the  progressions  of  the  melt  lines  of  the  Tempilaq 
coatings  as  the  result  of  aerodynamic  interference  heating.  The 
model  was  lightly  spray-coated  using  an  airbrush  and  allowed  to  dry. 
When  desired  flow  conditions  had  been  established  in  the  test 
section  of  the  tunnel,  the  movie  cameras  were  activated  as  the  model 
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was  injected  rapidly  into  the  flow  using  the  hydraulic  ram  facility 
of  the  NSWC/WOL  Hypersonic  Tunnel.  Details  of  the  tunnel  facility 
are  available  in  Reference  (12).  The  model  remained  in  the  flow  for 
test  times  of  typically  10-20  seconds  and  was  then  retracted.  Color 
movies  af  the  melt  line  progressions  were  taken  from  top  and  side 
views  using  two  Kodak  Cine  Special  TI  cameras  and  Kodak  16mm 
Ektachrome  Type  EF  7241  film.  The  setup  is  shown  schematically  in 
Figure  5.  Movies  were  taken  typically  at  24  frames  per  second. 

PHOTOGRAPHIC  COI^DITIONS 

A number  of  considerations  are  involved  in  obtaining  good 
photographic  data  in  this  paint  technique.  The  use  of  stroboscopic 
lighting  is  recommended  to  eliminate  extraneous  heating  of  the 
coatings  from  incandescent  lamps.  Radiant  heat  from  photo  flood- 
lights can  produce  considerable  error  in  the  determination  of  the 
heat-transfer  coefficient,  h.  This  error  may  be  especially  signifi- 
cant for  the  lower  rated  coatings  in  the  100-150°F  melting  temperature 
range.  In  these  tests  strobe  lights  were  mounted  flush  on  the  out- 
siue  of  the  tunnel  windows  to  cut  down  reflections  and  scattering  of 
the  incident  light  by  the  windows.  Fairly  uniform  lighting  of  the 
model  was  achieved.  There  was  no  illumination  other  than  that 
provided  by  these  strobe  lights,  which  were  synchronized  with  the 
framing  rates  of  the  cameras. 

EXPERIMENTAL  RESULTS 


DATA  REDUCTION 

Values  of  heat- transfer  coefficient,  h,  were  obtained  for  various 
fin-cone  configurations  by  the  somewhat  standard  data- reduction 
technique  for  this  method  (1) . The  16mm  color  movies  were  projected 
on  grid  paper  at  a convenient  geometric  scale  which  preserved  the 
clarity  of  the  melt  lines.  A stop-action  projector  v/as  used,  and 
successive  locations  of  melt-line  contours  were  drawn  on  the  grid 
model.  For  a given  contour  the  frame  number  was  recorded,  from  which 
a melt  time  was  deduced  based  on  the  camera  framing  rate  and  the  time 
of  initial  exposure  of  the  model  to  the  flow.  Errors  in  h due  to 
uncertainties  in  both  of  these  factors  are  discussed  in  a later 
section  of  this  report. 

ISOHEATING  CONTOURS 

Examples  of  the  reduced  heat- transfer  data  in  the  forms  of  lines 
of  constant  heat-transfer  coefficient,  so-called  isoheating  contours, 
are  shown  in  Figures  6 through  9.  Figures  6 and  7 show  side  and  top 
views,  respectively,  of  a flush-mounted  unswept  fin,  while  Figures 
8 and  9 show  similar  data  for  a flush-mounted  60°-swept  fin. 
Discussions  of  various  aspects  of  the  aerodynamic  interference 
heating  patterns  are  left  to  Reference  (4)  . This  report  concentrate." 
on  various  features  of  the  experimental  method  itself.  Inspection  of 


Baltakis 
Tunnel, " 


F.  P.,  "Performance  Capability  of  the  NOL  Hypersonic 
NOLTR  68-187,  Oct  1968 
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Figures  6-9  shows  one  advantage  of  the  method  to  be  that,  in  effect, 
the  entire  area  of  interest  is  "instrumented."  Thus,  a quantitative 
measure  of  heat  transfer  is  obtained  for  almost  all  of  the  inter- 
ference flow  field.  This  is  a recognized  advantage  of  thermal 
mapping  techniques,  in  general.  One  may  choose  to  utilize  these 
data  in  selecting  judicious  locations  for  discrete  sensors  in 
complementary  or  fo].low-up  tests.  The  phase-change  paint  method 
determines  both  the  severity  and  extent  of  the  interference  heating 
and  identifies  regions  of  extreme  gradients  that  might  hove  been 
missed  by  thermocouples  or  heat-flux  gages. 

UNCERTAINTY  IN  HEAT-TRANSFER  COEFFICIENT 


VALIDITY  OF  SEMI-INFINITE  SLAB  ASSUMPTION 

The  determined  value  of  heat-transfer  coefficient  is  subject  to 
uncertainty  due  to  uncertainties  in  the  numerous  input  parameters 
involved  in  both  the  experimental  and  analytical  parts  of  the  phase- 
change  paint  method.  Of  primary  importance  is  the  premise  that  the 
model  behave  like  a semi-infinite  slab  for  the  duration  of  a test 
run.  This  condition,  that  is,  that  the  test  time  be  less  than  the 
thermal  diffusion  time,  is  satisfied  if  Equation  (9)  is  satisfied. 

For  a Teflon  model  of  0.25-inch  wall  thickness  with  thermophysical 
properties  as  in  Table  1,  the  thermal  diffusion  time  is  about  80 
seconds.  On  areas  of  the  model  where  curvature  effects  come  into 
play.  Equation  (10)  is  applicable.  The  leading  edges  of  the  fins 
provide  the  most  severe  restriction,  for  which  the  thermal  diffusion 
time  is  about  20  seconds.  The  maximum  exposure  times  of  the  model 
in  the  flow  were  generally  about  20  seconds,  so  in  all  cases  the 
test  times  for  which  data  are  shown  were  less  than  the  thermal 
diffusion  time  for  the  Teflon  model.  The  two  embedded  thermocouples 
in  the  cone  ahead  of  the  fins  showed  no  temperature  rise  from  T. 

during  a run.  The  thermocouples  embedded  in  the  fins  exhibited 
gradual  temperature  rise,  especially  at  higher  Reynolds  numbers. 

This  indicated  that  sides  of  the  fins  did  not  behave  like  a semi- 
infinite slab  for  the  duration  of  a run.  Consequently,  the  fin  side- 
heating data  are  not  considered  to  be  as  reliable  for  two  higher 
Reynolds  numbers  included  in  Reference  (4) . 


SENSITIVITY  TO  ERROR  IN  THERMGPHYSICAL  PROPERTIES  AND  INITIAL  TIME 
Overall  uncertainty  in  h is  the  result  of  uncertainty  in  the 
thermophysical  properties  of  the  model  material,  uncertainty  in  time 
of  initial  exposure  to  the  flow,  and  uncertainty  in  the  temperatures 
^i'  *^aw  *^m*  I^scall  Equation  (8),  given  by 


from  which 


S = — /t 

/pcHc 
■— P .. 

v't 
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h 
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It  follows  that 


h 


1 

2 \p 


k 


t / 3 


(12) 


With  respect  to  the  first  three  terms  on  the  right-hand  side  of 
Equation  (12),  the  estimated  accuracies -3)  for  the  properties 
of  Teflon  are;  density  p,  one  percent;  thermal  conductivity  k, 

seven  percent;  and  specific  heat  C , five  percent. 

P 

To  investigate  the  error  produced  in  h by  uncertainty  in  time  of 
initial  exposure  to  the  flow  for  constant  properties,  consider 
Figure  10,  which  shows  Ah/h  versus  exposure  time,  t,  to  flow, 
with  uncertainty  in  time.  At,  a parameter.  For  increasing  exposure 
times  the  effect  of  uncertainty  in  initial  tim.e^  t^,  decreases.  For 

these  tests  At  was  estimated  to  be  about  0.2  second,  during  which 
time  the  model  passed  through  the  core  flow  boundary  layer.  Data 
based  on  model  exposure  times  greater  than  2 seconds  should  be  used 
in  order  to  minimize  At  effects,  as  indicated  in  Figure  10. 


SENSITIVITY  TO  ERRORS  IN  TEMPERATURE  INPUTS 

With  regard  to  the  A6/3  term  in  Equation  (12) , cne  method  of 
attack  is  to  estimate  uncertainties  in  the  three  temperatures  (T^, 

T , and  T ) and  then  to  calculate  A3/3.  Examples  of  such  an 
3W  m 

approach  may  be  found  in  Reference  (14) . This  method  is  not  pursued 
here.  Instead,  analytic  relations  have  been  worked  out  which  show 
what  relative  uncertainty  is  produced  in  h for  a given  uncertainty 
in  an  input  temperature.  In  each  expression  two  temperatures  are 
considered  known  and  constant,  while  the  third  is  allowed  to  vary. 
The  three  relationships  were  obtained  by  differentiating  Equation 
(6)  with  respect  to  the  varying  temperature.  (These  expressions  are 
not  known  to  have  appeared  in  the  literature.  An  attempt  at  one 
such  relationship  (Ref,  (15))  is  considered  incorrect  by  an  inverted 
term  and  an  incorrect  sign.)  The  three  equations  are  given  by: 


(13) 

Kaufman,  L.  G.  , II,  Ling,  J.. , and  Johnson,  A.  R. , "Exploratory 
Tests  Using  Temperature-Sensitive  Paints  to  Obtain  Hypersonic 
Heat  Transfer  Data  on  Spheres  and  on  Fin-Plate  Models,"  GAC 
RM-487,  Grumman  Aerospace  Corp. , Sep  1970 

(14) 

Woods,  J, , Jr.,  and  Carlson,  J.  H.,  "Application  of  Temperature 
Sensitive  Paints  for  Aerodynamic  Heating  Analysis,"  MDC  A1419, 
McDonnell  Douglas  Corp.,  Nov  1971 

(15) 

Edney,  B.  E. , Bramlette,  T.  T. , Ives,  J.,  Hains,  F.  D. , and 
Keyes,  J.  W. , "Theoretical  and  Experimental  Studies  of  Shock 
Interference  Heating,"  Report  No.  9500-920-195,  Bell  Aerospace 
Co. , Oct  1970 
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where  a common  factor  r(3)  is  defined  as 

1 


r(3) 


/? 

r(3)  E g 


e^  erfc  3 - ij 


(13) 


(14) 


(15) 


(16) 


SENSITIVITY  TO  ERROR  IN  ADIABATIC  WALL  TEMPERATURE 

Equation  (13)  shows  that  error  in  k due  to  error  in  T is 

proportional  to  T.  Thus,  error  in  h_may  be  minimized  by  choosing 
test  conditions  with  small  T,  where  T is  defined  in  Equation  (7) . 

In  addition,  test  conditions  with  small  (T./T^„)  help  minimize  error 

in  h.  This  small  (T./T  .)  condition  is  advantageous  in  all  three 

1 aw 

equations.  The  (1  - T./T  ) factor  in  the  denominator  was  purposely 

not  included  in  the  definition  of  F (3)  in  order  to  emphasize  the 
importance  of  low  (T^/T^^)  ratio.  Figure  11  shows  Equation  (13)  for 

three  different  values  of  T for  T./T  = 0.7.  These  values  are 

m X aw 

typical  of  conditions  for  the  results  shown  in  Figures  6-9  and  in 
Reference  (4) . 

SENSITIVITY  TO  ERROR  IN  INITIAL  MODEL  TEMPERATURE 

Equation  (14)  demonstrates  that  the  error  in  h varies  as 
(1  - T)  with  error  in  initial  temperature,  T^,  where  once  again  the 

initial  temperature  to  recovery  temperature  is  an  important  factor. 
Here  it  is  doubly  important  to  keep  the  ratio  (T./T  ) small,  the 

effect  being  a combined  one  of  (T./T  ) in  the  numerator  and 

X 3W 

[1  - (T./T  )]  in  the  denominator.  Figure  12  shows  Equation  (14) 

X dW 

for  three  different  values  of  T . 

m 


SENSITIVITY  TO  ERROR  IN  MELT  TEMPERATURE  OF  PHASE-CHANGE  COATING 

The  last  relationship.  Equation  (15) , specifies  what  uncertainty 
in  h results  from  uncertainty  in  melt  temperature,  T^^,  of  the  paint. 

?his  error  in  h is  proportional  to  the  ratio  (T^^/T^^)  and  again 
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varies  inversely  with  the  (1  - T./T^„)  factor.  Equation  (15)  is  shown 

X aW 

in  Figure  13  for  several  different  paint-melt  temperatures. 

GENERAL  COMMENTS 

The  analytic  expressions  in  Equation  (13)- (15)  should  provide 
guidelines  ror  the  design  of  experiments  using  the  phase-change  paint 
technique  in  order  to  minimize  undertainty  in  the  h values. 

Based  on  all  the  uncertainties  in  input  temperatures,  time,  and 
Teflon  properties,  most  of  the  h values  shown  in  Figures  6-9  appear 
to  be  good  within  about  20  percent. 

The  extreme  sensitivity  of  error  in  h due  to  uncertainties  in 
T^  or  in  Figures  12  and  13,  respectively,  reflect  the  fact  that 

the  solution  to  the  heat-conduction  equation  is  very  sensitive  to  its 
boundary  conditions,  which  themselves  predicate  1' ^ow ledge  of  T^  and  T^. 

Initial  temperature,  T^^,  is  generally  known  within  one  percent  from 

thermocouple  data,  and  coating  melt  temperature,  T^,  is  frequently 

known  within  one  percent  (see  Appendix  A) . Again,  the  uncertainty 
in  h goes  inversely  as  (1  - T./i\,  ) , so  a small  (T./T,  ) ratio  is 

X 3.W  X 3W 

always  desirable. 


CONCLUSIONS 

A marked  advantage  of  the  phase-change  paint  method  is  that,  in 
effect,  the  entire  area  of  interest  on  the  model  is  instrumented. 

In  addition,  models  generally  may  be  fabricated  faster  and  at  a 
lower  cost  than,  for  example,  thin-walled  models.  Even  if  flow 
conditions  in  a given  facility  do  not  allow  one  to  determine  h values 
with  low  relative  uncertainty  per  Equations  (13)-(15)  , it  .is  still 
possible  to  conduct  a relatively  inexpensive  initial  test  using 
paints  to  help  determine  placement  locations  for  discrete  sensors  for 
follow-up  tests. 

One  disadvantage  of  the  method  is  a problem  associated  with  most 
optical  methods,  namely,  the  tedious  reduction  of  the  photographic 
(movie)  data  to  visually  identify  the  melt-line  locations.  This 
identification  may  vary  somewhat  with  the  observer  and  with  the  clarity 
of  the  film  record.  In  addition,  camera  angle  and  model  attitude 
must  be  accounted  for  if  one  wishes  to  obtain  careful  spatial 
distributions  of  h over  the  model  surface  in  the  form  of  h versus 
distance. 

Overall,  the  phase-change  paint  technique  is  considered  capable 
of  providing  reliable  quantitative  heat-transfer  data.  The  technique 
is  especially  useful  on  complex  geometries  involving  aerodynamic 
interference  heating  patterns  of  initially  unknown  extent  and 
location.  With  careful  experiment  designing,  based  on  guidelines 
provided  here,  one  may  expect  to  obtain  heat- transfer  measurements 
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from  the  phase-change  paint  method  which  would  be  just  as  reliable 
as  results  from  other  methods  employing  discrete  sensors. 
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TABLE  1 

THERMOPHYSICAL  PROPERTIES 

3 

Density,  Ib/ft 

Specific  heat,  BTU/lb  °R 

Thermal  Conductivity, 
BTU/sec-ft-°R 

Melting  Temperature,  °F 

2 

Thermal  Diffusivity,  ft  /sac 


*See  References  (1) , (11)  and  (13) 
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FIG.  2 TEFLON  HEAT-TRANSFER  MODEL  (SCALE  IN  INCHES) 
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FIG.  3 SCHEMATIC  DIAGRAM  OF  FIN-COIME  MODEL 
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FIG.  4 SCHEMATIC  DIAGRAM  OF  FINNED  EXTENSION 


1.  KODAK  CINE  SPECIAL  II  CAMERAS  ; '02mm  LENS.  f2. 

2.  GENERAL  RADIO  STROBOTACS,  MODEL  1538  (HIGH 

3.  SELSYN  GENERATOR-DRIVEN  MOTORS 

4.  GENERATOR 

5.  CEC  CINE  SPECIAL  VARIABLE-SPEED  MOTOR 
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Mm  = 5 

Re  00 /FT  = 4.5  10^ 

GAP  = 0.0" 


FIG.  6 ISOHEATING  CONTOURS  FOR  FLUSH-MOUNTED 
UNSWEPT  FIN.  SIDE  VIEW 


23 


NSWC/WOL/TR  75-62 


M 00  = 5 

Re  00 /FT  = 4.5  x 10® 
GAP  = 0.0" 


h X lO^,  BTU/Ft2-SEC-0r 


FIG.  8 ISOHEATING  CONTOURS  FOR  FLUSH-MOUNTED  60°-SWEPT 
FIN.  SIDE  VIEW. 
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APPENDIX  A 

PEASE-CHANGE  PAINT  CALIBRATION 

A schematic  diagram  of  the  calibration  setup  is  shown  in 
Figure  A-1.  The  rig  consisted  of  a single  chromel-alumel  thermo- 
couple spot-welded  to  a 1/32-inch-thick  stainless  steel  plate.  The 
thermocouple  output  was  fed  through  an  electronic  cold  junction 
compensator  into  a strip  recorder.  The  strip  recorder  was  calibrated 
using  a millivolt  potentiometer  source.  The  plate  was  heated  by  a 
heat  lamp  on  the  side  opposite  the  thermocouple. 

A roughly  1/2-inch  by  1/4-inch  strip  of  the  paint  being  checked 
was  brushed  on  the  plate  about  1/4  inch  to  the  right  of  the  thermo- 
couple and  allowed  to  dry.  The  heat  lamp  was  turned  on  and  a 
temperature- time  history  was  recorded  for  the  thermocouple.  When  the 
paint  coating  was  observed  to  turn  clear  (melt),  the  scribe  on  the 
recorder  was  raised  and  released  to  mark  the  spot  on  the  history 
curve.  The  heating  rate  could  be  varied  somewhat  by  moving  the 
plate  relative  to  the  heat  lamp.  Jones  and  Hunt(l)  found  almost 
no  variation  of  melt  temperature  with  either  pressure  or  heating 
rate  for  Tempilaq  coatings.  The  calibration  results  for  a number 
of  paints  are  shown  in  Table  A-1  for  a heating  rate  of  about  4°F/sec. 
The  observed  melting  temperatures  of  the  coatings  agreed  well  with 
their  rated  temperatures.  The  manufacturer  specifies  Tempilaq 
paints  to  melt  within  il  percent  of  their  rated  output.  Of  the  nine 
coatings  tested,  five  were  within  one  percent,  one  was  within  two 
percent,  two  v;ere  within  three  percent  and  one  was  within  four 
percent. 
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FIG.  A-1  SCHEMATIC  DIAGRAM  OF  TEMPERATURE- SENSITIVE  PAINT  CALIBRATION  SETUP 
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APPENDIX  B 


GENERAL  SOLUTION  OF  HEAT  EQUATION 

The  one-dimensional  heat-conduction  equation  given  by  Equation 
(1) , subject  to  the  boundary  conditions  in  Equations  (2) , (3)  and 
(4),  has  a general  solution  (Refs.  (6),  (16))  for  temperature,  T,  at 
time,  t,  at  any  depth,  x,  into  the  semi-infinite  slab  given  by: 


where  a is  the  thermal  diffusivity  defined  as  a = k/pC 


P’ 


At  the  model  surface,  x = 0 and  T(0,t)  is  assumed  co  be  equal 
to  the  melt  temperature  of  the  phase-change  coating^,  so  that  T(0,t)  = 
Tj^.  The  left-hand  side  of  Equation  (B-1)  becomes  T as  defined  by 

Equation  (7) , and  the  right-hand  side  may  be  simplified  since 
erf(O)  = 0.  Equation  (B-1)  becomes 


T = 1 - 


(B-2) 


Recalling  the  definition  of  3 from  Equation  (8) , Equation  (B-2) 
becomes 

.2 

T = 1 - e*^  erfc  3 , 


in  agreement  with  Equation  (6) . 
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